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ABSTRACT 

Numerical solutions of mutually coupled time-dependent equations of con- 
tinuity , momentum and energy balance are presented to illustrate the effect of 
eddy turbulence on the neutral composition and temperature of the lower and 
upper thermosphere. The procedure adopted in this paper eliminates the nec- 
essity of making ad-hoc assumptions about the conditions at the turbopause level 
and allows the study of the two regions in a self-consistent manner. From the 
illustrative examples comprising parametric changes in the eddy diffusion co- 
efficient, the specific roles of eddy turbulence in the development of theoretical 


models of the thermosphere are discussed. 
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INTRODUCTION 

In developing theoretical models of the upper atmosphere, it is customary 
to divide the thermospheric region into two domains with a boundary in the alti- 
tude range of 100-120 km. This boundary is usually identified with the turbo- 
pause level which by definition is the altitude at which the eddy and molecular 
diffusion coefficients become equal. The two regions, therefore, are distin- 
guished according to the relative importance of molecular and eddy transport. 
Because of the differences in the geophysical conditions of the two regions, their 
studies are pursued independently and the link between the two regions is con- 
trived by making ad-hoc assumptions about the boundary values at the turbopause 
level. The boundary values thus chosen are not necessarily consistent with the 
equations of continuity , momentum and energy transport of this region and are , 
in fact, among the most uncertain elements in the development of theoretical 
models. 

An understanding of the nature of variation of neutral composition and tem- 
perature in the altitude region of the turbopause is needed not only for the gen- 
eral understanding of the thermospheric behavior but also for the understanding 
of a large number of ionospheric phenomena which are strongly influenced by 
changes in the neutral atmosphere. In a series of papers, Chandra and Herman 
[1969], Chandra and Stubbe [1971] and Chandra et al. [1972] have shown how the 
changes in neutral composition at the turbopause trigger a complex chain of 
events, which not only affect the distributions of both the ionized and neutral 
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constituents and their temperatures, but also the emission rates of OI 6300 A 
during a SARARC event. The concept of changing neutral composition has been 
introduced in several studies for explaining the seasonal and magnetic storm re- 
lated changes in the ionosphere [King, 1967; Duncan, 1969; Fatkullin, 1973]. 
Some of these variations in the lower thermosphere have also been inferred 
from OGO-6 neutral mass spectrometer measurements [Reber et al. , 1973). 

On a day-to-day basis , there is considerable variability in the conditions of the 
lower thermosphere. This has been inferred from several rocket measurements 
of density and temperature [Spencer et al. , 1970; Smith et al. , 1972]. Although 
no clear pattern about the nature of these variations has been established yet, it 
is obvious that the conditions in the lower thermosphere are not invariant as has 
been assumed in the development of several atmospheric models of the sixties 
[CIRA 1965, U.S. Standard Atmosphere 1966]. 

As a first step towards understanding the nature of variations in the altitude 
region of the turbopause, it is important that we treat the thermosphere as a 
single entity and solve the relevant equations of the two regions in a self- 
consistent manner. This will eliminate the necessity of using turbopause as a 
'lid' separating the two regions of the thermosphere. We must, in addition, 
understand the nature of eddy turbulence in the lower thermosphere which is a 
dominantly controlling factor for both the composition and thermal structure of 
this region. Unfortunately, no simple way is available at present to theoretically 
relate thermospheric turbulence in terms of basic geophysical parameters. The 
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effect of eddy turbulence is, therefore, studied by defining eddy diffusion and eddy 
conduction fluxes , analogous to molecular diffusion and molecular conduction fluxes 
of particles and energy, respectively [Colegrove et al. , 1965, 1966; Johnson and 
Wilkins, 1965; Johnson and Gottlieb, 1970]. The eddy diffusion coefficient 
which determines these quantities is estimated from the sodium vapor trail ex- 
periments [Blamont and de Jager, 1961; Zimmerman and Champion, 1963; 

Justus, 1967; Zimmerman et al. , 1971]. Its values have been measured over a 
wide range and it has not yet been possible to establish a clear pattern for its 
latitudinal, diurnal, seasonal and height variations . At present, therefore, the 
effect of eddy diffusion can be best studied by parameterizing the eddy diffusion 
coefficient, and estimating it indirectly by comparing the observational data with 
theoretical calculations. This procedure has been adopted by several workers 
[Ivel'skaya et al. , 1970; Shimazaki and Laird, 1970; Shimazaki, 1967, 1971; 
George et al. , 1972]. In most of these studies, however, the equations of 
continuity are solved with the temperature assumed as a prefixed or known pa- 
rameter. Since changes in eddy diffusion coefficient also affects the thermal 
structure significantly, it is clear that the effect of eddy mixing is much more 
complex than revealed by these calculations [Iwasaka, 1973]. 

The purpose of this paper is to discuss some of these complexities and their 
physical implications from the simultaneous solutions of the equations of con- 
tinuity and energy transport for the thermospheric region. These calculations 
lead to the development of self-consistent models of the neutral composition and 
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temperature of the lower and upper thermosphere and allow the study of specific 
roles of eddy turbulence in a wide variety of geophysical phenomena. 

BASIC EQUATIONS 


Equations of Continuity 

The neutral constituents considered in this study are O, 0 2 , and N 2 , the 
main constituents contributing to the energetics of the thermosphere in the alti- 
tude range of 85-450 km. In deriving their distributions, we have assumed N 2 
to be in diffusive equilibrium throughout the above altitude range. The distribu- 
tion of N 2 is, therefore, uniquely determined from temperature distribution 
which is obtained by solving the energy balance equation. The only motion at- 
tributed to N 2 is the 'breathing' motion caused by the expansion and contraction 
of the atmosphere. The distribution of N 2 and its velocity are thus given by the 
following expression: 


[N 2 ] = 


_ [NaL T o 


exp 


f Z — 

~L r n 2 


V N , = 


g 


f Z g(z') 9T (z») + H ± fp * 

J t^) at N2 at 


(i) 


( 2 ) 


where the suffix o corresponds to the value at the lower boundary z 0 (= 85 km). 
Here 


Pn 2 = partial pressure of N 2 = [N 2 ] kT (3) 

kT 

H N = N 2 scale height = — — - (4) 

2 ni N2 g 

where g is the acceleration due to gravity, k the Boltzmann constant, mn 2 the 
molecular weight of N 2 and T the temperature. 
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The assumption of diffusive equilibrium for O and O 2 is valid only in the 
upper thermosphere (z <; 120 km). In the lower thermosphere their distributions 
are controlled by chemical, photochemical and transport processes and must be 
determined by solving the appropriate equations of continuity as follows: 


= Qi - Lj - ( ni Vj ) , (i = 1, 2) (5) 

where qj and Lj are production and loss rates, Vj the vertical transport velocity 
and the suffixes i = 1 and 2 designate O and 0 2 , respectively. 

The terms qj and Lj arise from photodissociation of 0 2 [in the Schumann- 
Runge Continuum (X = 1250-1700 A) and Schumann-Runge Band (X = 1700- 
2100 A)] and from 3-body recombination and are calculated here in the same 
way as discussed by Shimazaki [1971]. The velocity Vj in eqn. (5) consists of 
3 terms: (i) the velocity V N ^ of the major constituent N 2 , (ii) molecular diffu- 
sion velocity V jD , and (iii) the eddy diffusion velocity V iK , relative to N 2 . Thus 

= V N2 + V iD + V iK (6. a) 

The vertical mass motion of the neutral gas is thus given by the mean 
value : 


£ ^ Vi 3 

v z = — , Pi = m, ni , p = £ Pj 

P i=l 


(6.b) 


where the suffixes i = 1, 2, and 3 denote O, 0 2 and N 2 , respectively. Here 
V N2 is given by equation (2) and Vj D and Vj K are given by the following 
expressions: 
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V iD = -D 


i3 


V iK = - K 


' 1 

9n, 

i 

+ T 

3T 

. n i 

9z 

0 Z 

" 1 

9ni 

+ 1 

0T 

. n i 

9z 

T 

9z 


+ 


Hj 


H 


(7. a) 
(7.b) 


where Dj 3 and K are molecular and eddy diffusion coefficients relative to N 2 , 
respectively, Hj is the respective scale height and H is the average scale height 
which approximately corresponds to the scale height of N 2 . The coefficient Dj 3 
is calculated from the expression of a simple binary diffusion coefficient as 
follows [Chapman and Cowling, 1952]: 


Di3 = 


kT • mj + m N2 | 


, (i = 1, 2) 


( 8 ) 


8 °i3 2 { ni + [ N 2]} L 27rm * mN 2 
where o i3 is the collision diameter ( a 13 = 2.4x 10' 8 cm, o 23 = 3.4 x 10 _8 cm) 

[Lettau, 1951]. 

An expression for K in terms of simple geophysical parameters is not 
available, as mentioned earlier. Its functional form, therefore, must be de- 
scribed in terms of assumed parameters. In this paper, we ha.ve adopted a 
model of K(z) given by Shimazaki [1971] , which basically relates the observa- 
tional data in terms of a few simple parameters . Denoting the maximum value 
of eddy coefficient by K M , and the corresponding height by z m , we write: 


K(z) =* 
where S 


K m exp [ -S j (z - z m ) 2 ] , z >z m 

/ 

1<K M " K 0 ) exp [-S 2 (z - z m ) 2 ] +K 0 exp [s 3 (z - z m )] 
i, S 2 and S 3 are shape parameters and K 0 is a constant. 


( 9 ) 

> Z 5: z m 
With suitable 



choice of these parameters, the eddy diffusion coefficient can be made to fall 
off rapidly above the height of the maximum, simulating the condition of rapid 
cessation of turbulence above the turbopause. 

Energy Balance Equation 

The energy equation for determination of thermal structure above 120 km 
has been discussed in detail by Chandra and Sinha [1973] . In this paper this 
equation has been modified to include the terms arising from eddy transport 
and 3-body recombination. The basic form of energy balance equation solved 
in this paper is as follows: 

Z] (Wj p j " itO = Qi + Qd + Qr + Qc + Q e ~ Qe ( 10 ) 

j = l ' ' 

d a „ 3 

with dt = si + Vz si ' 

j = 1, 2, 3 for O, 0 2 , N 2 , respectively, and Cpj is the corresponding specific 
heat at constant pressure. 

The terms on R.H.S. of eqn. (10) represent the contributions from various 
sources and sinks defined as follows: 

Q i = Thermal energy from photoelectrons produced by photo-ionization of 
O, 0 2 and N 2 in the EUV range of solar spectrum. 

Q d = Energy available from photodissociation of 0 2 in the Schumann-Runge 
continuum and Schumann-Runge band. 

Q g = Heat loss arising due to infrared radiation (63ju) emitted by thermally 


excited atomic oxygen. 
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The expressions for Qj , Q d and Qg are given in a previous paper [Chandra 
and Sinha, 1973] and are calculated here in exactly the same way. Q r is the 
energy input rate due to 3-body recombination of O: 

O + O + M -* 0 2 * + M (11. a) 

0 2 * 0 2 + hi> H (11. b) 

The recombination of O normally results in the production of 0 2 in the ex- 
cited state (denoted by the superscript *) and its subsequent de-excitation 
through emission in the Herzberg band (frequency i> H ). The recombination 
energy is calculated from the following expression: 

Q r =«N [O] 2 (e d - e H ) (12. a) 

where e d is the dissociation energy of 0 2 (=5. lev), e H the average energy of the 
Herzberg band (=4.25ev), N the total number density and a the rate coefficient. 
a. has a strong temperature dependence and is given by [Campbell and Thrush, 


1967]. 


a = 3 x 10 


-33 


300 


- 2.9 


cm 6 sec 1 


(12. b) 


Finally, Q c and Q e represent contributions due to thermal (molecular) and 
eddy conduction, respectively: 


Qc = ai 


, 0T 
X 


(13. a) 


Qe - 


pkcp (f + r 


(13. b) 


where X is the thermal (molecular) conductivity ,p the total density, Cp the 


specific heat at constant pressure, and T is the adiabatic lapse rate ( T = 
of the neutral gas . 


JL 

Cp 
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Numerical Procedure and Boundary Conditions 

The equations of continuity and energy balance discussed in the preceding 
sections are formally all of the same type . They are partial differential equa- 
tions of first order in t (time) and second order in z (height). The numerical 
solutions of these equations are obtained by using an implicit finite difference 
technique, with integration steps At = 30 minutes and Az = 1km. The scheme 
for solving these equations is basically the same as discussed by Herman and 
Chandra [1969] and Stubbe [1970]. 


Equation (5) is relevant only in the lower thermosphere where O and 0 2 are 
minor constituents. In the upper thermosphere O becomes a major constituent 
as a result of diffusive separation and the distribution of individual constituents 
and their velocities can be calculated from appropriate equations analogous to 
(1) and (2), respectively. The upper boundary for equation (5), therefore, has 
to be chosen at a height where N 2 is still a major constituent and above which 
level molecular diffusion dominates over eddy diffusion, production and loss 
terms. We may define this height as a transition height z t , to distinguish it 
from a turbopause height. 

The boundary conditions at z t may then be determined from the requirement 
of the continuity of flux under the specified simplification (K < D i3 ): 


3ni 1 cIT _1 " 

nj dz T dz + Hj 



H -Ei 
Hl 3t 


(i = 1, 2) 
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If z t is chosen above 120 km, where the typical values of D i3 ^ 10 8 cm 2 sec -1 , 
the expression on the R.H.S. becomes negligibly small. In this case, the 
boundary values at z t may be simply determined from the following condition 
representative of diffusive equilibrium distribution: 


On, 

9z 


z = z. 


1 9T 1 
n i l T 9z + H ; 


(14. a) 


z = z. 


“t " t 

For the lower boundary z Q , we may assume the photochemical terms to 
dominate over the transport term. The time-dependent boundary conditions for 
O and 0 2 densities in this case may therefore be determined from the following 
equation: 


9nj 

9F 


= (Qi - L i) 


z = z„ 


i.e. , 


n j (z 0 , t) = n; (Z 0 , t - At) + At (qj - Lj) (14. b) 

z o > 1 

The concept of the transition height z ( is useful only for solving the contin- 
uity equations of O and 0 2 • For solving energy balance equation (eqn. 10), it is 
not necessary to treat the lower and upper regions of the thermosphere differen- 
tially through introduction of a transition height . The upper and lower boundary 
conditions used in solving equation (10) are as follows: 


At 

and at 


z = z u (450 km): 
z = z Q (85km): 


9T 

9z 


T(z 0 ) 


= 0 

= 180°K 


r ( 14 - c ) 
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NUMERICAL RESULTS AND DISCUSSION 

In presenting the numerical solutions of the energy equation, the choice of 
solar flux is clearly of critical importance. In a recent paper, Hinteregger 
[1970] has given a table of spectral distribution of solar EUV flux which he be- 
lieves to be representative of medium solar activity. The values of EUV flux 
given in this paper differ considerably from those given by Hinteregger et al. 
[1965] for reportedly low solar activity and are, on the average, a factor of two 
less than their earlier values. Hinteregger [1970] also believes that the solar 
flux in the Schumann-Runge Continuum as tabulated by Hinteregger et al. [1965] 
is much too high and should be reduced by a factor of 3. Using the solar flux 
data given by Hinteregger [1970 ] and the currently accepted range of values of 
eddy diffusion coefficient we find the exospheric temperature is much too low for 
the conditions of medium solar activity (T^ < 600° K). This is in agreement 
with the conclusion arrived at by Roble and Dickinson [1973] . On the other 
hand, the use of the solar flux values given by Hinteregger et al. [1965] gives an 
exospheric temperature much too high for low solar activity (Too >1500°K) con- 
sistent with our earlier result [Chandra and Sinha, 1973] . Because of the para- 
metric nature of our present study, we have adopted solar flux values which can 
accommodate large variations in eddy diffusion coefficient and yet give the tem- 
perature variations in a reasonable range . The flux values in EUV range used 
in this paper are the same as given by Hinteregger [1970]. The flux values in 
the Schumann-Runge Continuum are taken from Hinteregger et al. [1965] and 
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the flux values in the Schumann-Runge Band are the same as used by Shimazaki 
[1971]. The absorption and ionization cross sections are appropriately taken 
from Hinteregger et al. [1965] and Shimazaki [1971]. 

The eddy coefficient (K-model) as defined by eqn. (9) has several param- 
eters which affect the altitude distribution of the eddy diffusion in the lower 
thermosphere. For the purpose of parametric studies, however, here we con- 
sider only the variations in K M and z m , which are the main parameters describ- 
ing the general features of the eddy turbulence in this region. The overall 
changes in K(z) associated with the changes in K M and z m are shown in Figures 
(la) and (lb), respectively. In both these plots, the values of Sj, S 2 , S 3 and 
K 0 as defined in eqn. (9) are fixed as follows: 

Sj = 0.05km" 2 , S 2 = 0.03km" 2 , S 3 = 0.07km" 1 , Kq = 10 6 cm 2 sec _1 . 

For the case depicted in Figure (la), z m has a fixed value of 110km and 
K m varies from 10 6 to 10 8 cm 2 sec" 1 . The effect of this change on K(z) is most 
noticeable in the immediate vicinity of the height of the maximum z m . The var- 
ious profiles merge into a single profile about 15 km below z m . In the altitude 
region above z m , K diminishes to a low value well before the transition height 
z t (=130 km) for all the four K M values shown. Clearly here higher K M value 
corresponds to higher K-value at all altitudes above 95 km. In contrast to this 
situation the variation of z m from 120 to 105 km with K M fixed at 5 x 10 7 cm 2 
sec' 1 (Fig. lb), affects the K-values differently in different altitude regions. 
The lowering of z m yields a higher K-value in the lower altitude region and 
lower K-values in the higher altitude region. 
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As a direct consequence of the altitude variation in K(z), the eddy term as- 
sociated with the adiabatic lapse rate in equation (13. b) acts as a heat source or 
sink in the various altitude regions. Denoting this term by Q ad we may write 
from eqn. (13. b). 


Qad 


gpK 


' 1 3K 1 9p ‘ 
K dz + p 9z 


(15) 


Since 


1 dp 
p dz 


1 , 13K . n . 

‘ H> and K 31 < ° forz>z m’ 


Q ad is always negative in the altitude region above the height z m and hence is an 


effective heat sink in this region. In the altitude region below z m , 17 > 0 and 


_1 9K 19 P 
[K 9z + p 3 zJ* 


Q ad may act either as a heat source or sink depending upon the sign of 

The variations in K(z), therefore, results in the increase or 
decrease of the neutral temperature in the various regions of the thermosphere 
depending on the nature of Q ad . The eddy related changes in the neutral compo- 
sition are even more difficult to describe in general, since the changes in K(z) 
affect the effective scale height of the minor constituents both directly (via 
mixing) and indirectly (via neutral temperature) . Some of the intricacies of 
these interacting processes can be described in detail only through time- 
dependent numerical solutions presented in the following sections. 

Variation of K M (z m Constant) 

The temperature and atomic oxygen density profiles corresponding to the 


K-models discussed in Figure (la) are shown in Figure (2) for the local noon 
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condition. The numerical results are shown here only up to about 200km since 
there are no basic changes in the trends of the various profiles above this alti- 
tude region. The temperature in all the cases increases continuously with alti- 
tude until it becomes nearly isothermal (at about 300 km), and the atomic oxygen 
density, in accordance with the diffusive equilibrium distribution, falls off ex- 
ponentially at a rate determined by its scale height. Both the temperature and 
atomic oxygen profiles have several interesting features in the altitude region 
of 85-200 km, which are summarized as follows: 

(1) The temperature varies inversely as K(z) according to the qualitative 
deductions discussed in the preceding section. 

(2) The temperature, in general, does not increase monotonically with al- 
titude , but instead goes through a maximum and a minimum in the alti- 
tude region of 95-120 km. The change from maximum to minimum de- 
pends upon the values of K(z) and is attributed to the heating and cooling 
associated with the adiabatic lapse rate given by eqn. (15). Low K M 
value, however, yields a monotonic temperature profile, as exemplified 
by K m = 10 6 cm 2 sec _1 curve. 

(3) The variations in K M do not have any significant effect on [O] in the 
altitude region below 100 km, since this region is basically under 
photochemical control. In the altitude region above 100km, [O] de- 
creases with increase in K(z). The effect of increasing K(z) is to de- 
crease the effective scale height of atomic oxygen directly (through in- 
creased mixing) and indirectly (through a reduced neutral temperature). 
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The diurnal variation of the exospheric temperatures (taken at 450 km) for 
representative values of K M in this series is represented in Fig. (3). As seen 
from this plot, the change in K M affects both the diurnal amplitude and the ab- 
solute values of temperature but not the basic diurnal characteristic with respect 
to the times of the maximum and minimum. A decrease in K M from 10 8 to 10 6 
cm 2 -sec -1 results in increase of the exospheric temperature by about 50%. A 
change of this magnitude is comparable to the changes associated with the solar 
activity over a solar cycle. 

In presenting the numerical results of Figs . (2) and (3) the eddy diffusion 
models were assumed to be independent of time, since their diurnal variations 
are not known with any certainty. In the spirit of parametric studies, we have 
investigated the effect of varying K M over a diurnal time scale. The time vary- 
ing models of K M have little effect on the altitude and diurnal characteristics of 
temperature and densities, though they do change their absolute values. A de- 
crease in K m by an order of magnitude from day to night , for example , in- 
creases the exospheric temperature by about 20% without changing its basic 


altitude characteristics. 
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Variations of z m (K M Constant) 

The altitude variations in temperatures and [O] as a result of changing z m 
are quite different from those due to the variations in K M . Fig. (4) is a typical 
example of results for changing z m from 120 to 105 km with K M = 5 x 10 7 cm 2 
sec' 1 (Fig. lb) during the day with one order of magnitude reduction in K-value 
during the night. (As mentioned already, the results for K M = 5 x 10 7 cm 2 sec _1 
throughout the day and the night are basically similar except for some differ- 
ences in the absolute magnitudes.) Unlike the changes associated with K M , the 
changes in z m have different effects in the upper and the lower thermospheric 
regions. The reduction in z m results in the increase of the temperature in the 
higher altitude region and decrease in the lower altitude as a direct consequence 
of the redistribution of K(z) and hence of the adiabatic lapse rate related cooling 
(heating) rate (eqn. 15). 

The increase in K(z) in the lower thermosphere reduces the atomic oxygen 
in this region. However, this reduction does not persist at higher altitudes be- 
cause of the compensating effects from the increased temperature. The cross- 
over in temperature from a lower value in the lower thermosphere to a higher 
value in the upper thermosphere produces a similar effect on atomic oxygen 
and other neutral constituents through corresponding changes in their scale 
heights. More specifically, by lowering of z m from 120 to 105km, the atomic 
oxygen density decreases in the lower thermosphere due to reduction in its scale 
height. This decrease persists up to about 200-350 km, above which the pattern 
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gradually changes to increase in [O] . This kind of cross-over occurs in other 
constituents also. The altitude at which the crossover occurs for a specific 
change in z m depends on the constituent. For the heavier constituents, (0 2 , N 2 ) 
and also for the total density (Fig. 5), the crossover takes place in a much lower 
altitude range. An example of the changes in 0/N 2 and 0/0 2 associated with the 
change in z m is also shown in Fig. (5). The effect of decreasing z m , as seen in 
this plot, is to decrease 0/0 2 and 0/N 2 . This decrease is much more apparent 
in the higher altitude region (z > 160 km) and is a direct consequence of the dif- 
ferences in the scale heights of the three constituents above the transition height. 
SUMMARY AND CONCLUDING REMARKS 

In this paper, we have attempted to emphasize the importance of coupling 
between the lower and the upper thermosphere and the specific effects of eddy 
turbulence on the thermal structure and the neutral composition in the thermo- 
spheric region. Because of the complex nature of interactions of the various 
processes, we numerically solved the time-dependent equations of continuity 
and energy balance of the relevant neutral constituents in a self-consistent man- 
ner and without making ad-hoc assumptions about the composition and tempera- 
ture changes at the turbopause. 

Our study relating to the effects of eddy turbulence covered a wide range of 
variations in the eddy diffusion coefficient without relating these variations to any 
specific geophysical causes. The numerical results presented in this paper cor- 
respond to the geographic equator, though the qualitative inferences drawn from 



20 


these calculations are quite general and are applicable to the high and mid- 
latitude regions. Notwithstanding the parametric nature of this study, it seems 
reasonable to draw certain conclusions regarding the geophysical implications of 
the types of results obtained. Variation in the magnitude K M and altitude z m of 
eddy coefficient may have important implications in regard to seasonal variation 
of the neutral temperature and composition [Kockarts, 1972]. 

In particular, the temperature and composition changes associated with 
variations in z m are analogous to those actually observed during a geomagnetic 
storm at middle and high latitudes. A large reduction in 0/N 2 ratio, together 
with increase in exospheric temperature, are among the characteristic changes 
in the thermosphere during a magnetic storm. [Jacchia et al. , 1966; Blamont 
and Luton, 1972; Reber et al. , 1973] . We observe that lowering of induces 
precisely these types of effects. Thus we are led to the tentative conclusion 
that during the main phase of a geomagnetic storm the turbulence in the lower 
thermosphere is shifted to lower altitudes . This is physically plausible in view 
of low altitude deposition of energy in the auroral region due to energetic particle 
precipitation and Joule dissipation during magnetic storms. Subsequent trans- 
port of energy to the mid-latitude regions is expected to yield enhanced turbu- 
lence at lower altitudes, which is schematically analogous to lowering of the 
z m -value. 

Obviously, by the very nature of the process, eddy turbulence is generally 
a highly random and variable phenomenon. Variations in the eddy diffusion 
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coefficient, a purely phenomenological parameter, are to be expected with 
changes in the geophysical conditions, such as tropospheric and stratospheric 
disturbances of meteorological nature, gravity and tidal waves, local solar 
heating , magnetospheric and field-aligned perturbations , etc . The effect of 
such variations on thermospheric temperature and composition can, in principle, 
be incorporated in model calculations via appropriate modifications in the eddy 
coefficient, as illustrated here. It is, therefore, evident that inclusion of the 
eddy turbulence as a basic physical process is important for constructing a 
coherent model of the upper atmosphere and the ionosphere. 
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Figure 1 
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Figure 5. 


FIGURE CAPTIONS 

a. The altitude profiles of K(z) for K M varying from 10 6 to 10 8 
cm 2 -sec -1 . 

■b. The altitude profiles of K(z) for z m varying from 105 to 120km. 
The altitude profiles of temperature and atomic oxygen density 
for K-models shown in Figure 1-a. 

The diurnal variations in exospheric temperature for varying 
values of K M . 

The altitude profiles of temperature and atomic oxygen density 
for K-models shown in Figure 1-b. 

The altitude profiles of 0/N 2 , 0/0 2 and the total density (mass) 
for K-models of Figure 1-b. 
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